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ABSTRACT 


The  intermittent-eddy-viscosity  hypothesis  is  proposed  for  the 
outer  layer  in  a  turbulent  boundary  layer  flow,  and  it  assumes  an  alter¬ 
nate  appearance  of  zero  eddy  viscosity  in  the  non-turbulent  region  and 
of  a  constant  eddy  viscosity  in  the  turbulent  region.  The  equation  of 
motion  is  simplified  by  assuming  the  velocity  defect  law  and  that  the 
velocity  defect  is  small  compared  with  the  free  stream  velocity.  With 
these  assumptions  the  distributions  of  mean  velocity  and  shear  stress 
in  the  outer  layer  of  a  turbulent  boundary  layer  on  a  flat  plate  for  con- 
&ant  pressure  are  obtained  using  the  intermittency  factor  measured  by 
Klebanoff.  The  calculated  results  show  good  agreement  with  those 


measured. 


INTRODUCTION 


A  turbulent  boundary  layer  is  divided  into  two  parts ,  that  is ,  the 
wall  region  and  the  outer  region,  the  boundary  between  them  being  at 
about  fifteen  per  cent  of  the  boundary  layer  thickness  from  the  surface. 

It  is  a  physical  characteristic  of  the  wall  region  that  the  shear  stress 
is  almost  constant.  The  flow  adjacent  to  the  wall  is,  however,  laminar, 
whereas  the  flow  is  entirely  turbulent  in  the  layer  further  than  yUt/0=3O, 
where  the  mean  velocity  distribution  is  well  described  with  the  so-called 
logarithmic  wall  law.  Thus  the  wall  region  or  the  constant  shear  layer 
is  subdivided  into  three  layers,  the  laminar  sublayer,  the  transition 
region  and  the  turbulent  region  governed  by  the  logarithmic  law. 

The  physical  features  in  the  outer  region  are  that  the  mean  velocity 
distribution  has  similarity  if  it  is  plotted  in  the  form  of  velocity  defect, 
which  is  the  difference  from  the  free  stream  velocity,  instead  of  velocity 
itself.  The  velocity  defect  law  in  the  logarithmic  form  can  describe  the 
velocity  profile  well  into  the  wall  region,  but  it  can  describe  the  velocity 
in  only  a  small  part  of  the  outer  region.  The  power  law  for  the  velocity 
profile  is  purely  empirical,  though  it  describes  the  profile  fairly  well  for 
its  simplicity.  The  parabolic  distribution  proposed  by  Hama1  describes 
the  profile  all  over  the  outer  region.  The  wake  function  proposed  by 

O 

Coles  was  derived  empirically  as  an  expression  for  the  deviation  of 
the  velocity  from  the  logarithmic  wall  law.  As  the  velocity  distribution 
by  Coles,  however,  does  not  go  smoothly  into  the  free  stream  velocity, 
an  outer  wall  function  has  been  proposed  by  Cornish^  to  improve  this  de¬ 
fect  of  the  Law  of  Wake .  These  laws  are  useful  for  practical  purposes, 
although  they  seem  to  lack  physical  foundations. 

Assuming  that  the  eddy  viscosity  is  constant  in  the  outer  region, 
the  velocity  profile  in  the  form  of  a  co-error  function  was  obtained 
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analytically  by  Townsend4  and  by  Hinze^.  The  profile  shows  appreciable 
deviations  from  the  measured  values  at  the  outer  edge  of  the  boundary 
layer.  The  reason  for  the  discrepancy  may  be  due  to  the  characteristics 
of  flow  changing  intermittently  from  turbulent  to  non-turbulent  flow  in 
the  outermost  layer.  A  velocity  profile  was  derived  by  the  writer®  taking 
the  intermittency  into  account,  but  it  too  did  not  agree  well  with  the 
measured  results . 

The  outer  region  ranges  over  about  eighty-five  per  cent  of  the  turbu¬ 
lent  boundary  layer,  and  it  is  not  affected  directly  by  the  shear  at  the 
wall,  whereas  the  mean  flow  in  the  wall  region  is  governed  by  the  wall 
stress.  It  is  a  characteristic  of  the  layer  that  the  eddy  viscosity  is 
nearly  constant  within  a  portion  about  forty-five  per  cent  of  the  total 
thickness  from  the  wall  and  the  intermittency  is  observed  outside  of  this 
region. 

In  the  present  paper  the  mean  velocity  distribution  in  the  outer 
region  is  derived  analytically  taking  the  characteristics  of  the  flow  into 
account  and  is  compared  with  experimental  results . 

In  Chapter  1 ,  the  Intermittent-Mean-Velocity  Hypothesis  is  pro¬ 
posed,  and  in  connection  with  this  treatment,  a  new  mean  velocity  dis¬ 
tribution  near  the  center  of  a  turbulent  pipe  flow  is  derived  by  making 
use  of  the  mixing  length  according  to  Kikuradse' . 

In  Chapter  2 ,  the  Intermittent-Eddy-Viscosity  Hypothesis  is  pro¬ 
posed,  and  the  mean  velocity  distribution  is  derived  in  the  form  of  velocity 

Q 

defect  by  making  use  of  the  intermittency  factor  according  to  Klebanoff  . 

The  wake  function  by  Coles  and  the  outer  wall  function  by  Cornish 
are  compared  with  the  results  calculated  by  the  present  forumla  for  velocity 
distribution.  The  shape  factor,  eddy  viscosity  and  mixing  length  are  dis¬ 


cussed. 
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NOTATION 

Constants  in  logarithmic  velocity  defect  law  (Equation  2-27) 

Constants  in  the  expression  of  mixing  length  in  pipe  flow 
(Equation  1-10) 

Constant  in  the  logarithmic  wall  law  (Equation  1-14) 

(UB S/£g ) (d $/dz)'  constant  in  Equation  2-^13 
Shape  factor 
Mixing  length 
Pipe  radius 

Eddy  Reynolds  number  based  on  £.  and  respectively 

Mean  velocity  in  x-direction 

Free  stream  velocity  or  velocity  at  pipe  center 

Mean  velocity  in  the  turbulent  region  of  intermittent,  outer  layer 

Friction  velocity,  Ut  =  . 

Mean  velocity  in  y-direction 

Turbulent  velocity  fluctuations  in  x-  and  y-directions ,  respectively 

Turbulent  shear  stress  per  unit  mass 

Distance  measured  along  surface 

Distance  normal  to  surface  measured  from  surface 

Constant  in  the  expression  of  eddy  viscosity  by  Clauser 

Intermittency  factor 

Boundary  layer  thickness 

Displacement  thickness 

Mean  eddy  viscosity  in  outer  layer,  £  —  f£. 

Eddy  viscosity  in  the  turbulent  region  of  outer  layer 

S  —  SW$ -o.78) 

Momentum  thickness 

Constant  in  the  logarithmic  wall  law,  Equation  1-14 


\  Constant  in  the  expression  of  the  center  velocity  in  pipe  flow 
if  Kinematic  viscosity 

5  5  -  K/S 

j>  Fluid  density 

T  Shear  stress 

X0  Shear  stress  at  wall 


CHAPTER  1 


THE  INTERMITTENT- MEAN -VELOCITY  HYPOTHESIS 
1-1.  Introduction 

It  was  shown  by  Schubauer®  that  the  turbulent  energy  distribution 
in  a  boundary  layer  divided  by  the  intermittency  factor  coincides  closely 
with  the  distribution  in  pipe  flow.  According  to  Schubauer,  such  close 
agreement  could  not  occur  unless  the  average  turbulent  energy  was  dis¬ 
tributed  through  the  turbulent  regions  of  a  boundary  layer  just  as  it  is 
in  the  fully  turbulent  section  of  a  pipe  flow. 

This  fact  suggests  that  the  mean  velocity  distribution  in  the  turb¬ 
ulent  region  in  a  boundary  layer  may  coincide  with  that  in  a  fully  turbu¬ 
lent  pipe  flow,  because  the  turbulent  energy  is  taken  from  the  mean-flow 
kinetic  energy  through  the  turbulent  shear  stress  in  that  layer.  From 
careful  study  of  the  oscilloscope  records  of  velocity  fluctuation,  it  was 
noticed  by  Klebanoff8  that  "the  trace  has  somewhat  of  a  square-wave 
appearance  in  the  intermittent  region,  and  that  the  non-turbulent  regions 
seem  to  be  at  a  constant  level  corresponding  to  that  of  the  free  stream, 
while  the  turbulent  regions  are  seen  to  be  centered  about  some  lower 
level.  The  difference  between  the  velocity  of  the  outside  potential 
flow  and  that  existing  in  the  turbulent  regions  seemed  to  depend  on  how 
far  past  the  measuring  position  the  instantaneous  edge  of  the  layer  ex¬ 
tended  at  the  particular  instant". 

If  the  mean  velocity  level  in  the  turbulent  region  may  be  assumed 
to  be  expressed  with  the  extrapolation  to  the  outer  layer  of  the  logarithmic 
wall  law  in  the  deeper,  non-intermittent  turbulent  region,  which  is  valid 
in  the  fully  turbulent  pipe  flow,  the  average  mean-velocity  distribution 
in  the  intermittently  turbulent,  outer  layer  will  be  described  by  the  average 
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of  the  free  stream  velocity  and  that  of  the  logarithmic  wall  law,  that  is. 


u  —  (/-r) u.+  rut  , 


(i-0 


where  U,  Uo  and  Ut  denote  the  average  mean  velocity,  the 
free  stream  velocity  and  the  mean  velocity  in  the  turbulent  region  accord¬ 
ing  to  the  logarithmic  law,  respectively,  and  <f  denotes  the  intermittency 
factor.  The  expression  can  be  written  alternately  in  the  form  of  velocity 
defect,  which  is  known  to  be  appropriate  to  describe  the  mean  velocity 
distribution  in  the  outer  region, 


u,-u  —  r(u.-L4) , 


U-Z) 


The  mean  velocity  distribution  was  calculated  by  the  writer6  making  use 
of  the  following  empirical  formulas  by  Klebanoff  for  the  velocity  defect 


and  for  the  intermittency  factor, 


LkZ_Ut, - *7  sltf,£+2.3 

f  =4r(l-erfS )  -  1 


Uc-U 

Ur 


(1-3) 


5=s(i-o.7t),j  ('~4) 

7^  ==  ’  (IS) 


where  Ut  denotes  the  friction  velocity  1 [%/f  }  T0  the  wall  shear 
stress,  £  ,  density,  ,  the  distance  from  the  wall,  and  $  ,  the 

boundary  layer  thickness.  The  velocity  distribution  calculated  by 
Equation  (1-5)  shows  poor  agreement  with  the  experimental  results  as 
shown  in  Fig.  1. 
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Now  it  is  recalled  that  the  turbulent  intensity  in  a  boundary 
layer  divided  by  coincides  with  that  in  pipe  flow.  Then  the  equa¬ 
tion  for  (U0  —  Uf)/Ux  should  be  presumably  that  for  pipe  flow 
instead  of  for  boundary  layer  flow,  that  is 


Uc-Ut 

Uv 


-$.(>2  lofto-%-  + 


U-6) 


which  was  derived  by  Hinze^  based  on  the  measurements  by  Laufer^. 

In  this  case,  Uo  denotes  the  maximum  velocity  at  the  center  of  pipe. 

The  use  of  this  formula  does  not  give  good  results  either.  It  has  been 
noticed,  however,  that  for  y/S  >0. 15  the  logarithmic  velocity  dis¬ 
tribution  deviates  from  the  actual  velocity  distribution.  A  correction 

12 

function  for  this  deviation  has  been  empirically  introduced  by  Millikan 1  , 

13  3 

Hinze  ,  and  Cornish  . 

In  the  present  chapter,  a  new  formula  is  derived  analytically 
for  the  mean  velocity  distribution  in  a  fully  turbulent  pipe  flow,  and 
the  mean  velocity  distribution  in  a  boundary  layer  will  be  calculated 
using  the  new  formula. 


1-2  .  The  mean  velocity  distribution  in  a  fully  turbulent  pipe  flow 

According  to  the  momentum  transfer  theory,  the  turbulent  shear 


stress  can  be  expressed  as 


r-rff f-f- 


0-7) 


where  T  and  t  denote  shearing  stress  and  mixing  length,  respec¬ 
tively.  As  the  shearing  stress  in  a  pipe  flow  is  distributed  linearly, 


it  is  written  as 

r  =  T0  (I—  y/R) 


d-s) 


f 
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where  ft  is  the  radius  of  a  pipe.  According  to  Nikuradse7,  the  mixing 
length  in  a  pipe  flow  is  represented  by  the  expression 

--  —  o,  /4  -  O.oS  (I-  -fcf-  o.o6(l~  (j~f) 

Now  the  mixing  length  is  assumed  as 

1=  a 2  -  42(i--$-f- cz(l- 

where  the  numerical  constants,  a,  ,4",  and  c  are  left  to  be  deter¬ 
mined  later. 

(i)  Near  the  wall,  or  y/R  « l 

In  this  case,  the  mixing  length  is  known  to  be  proportional  to 
the  distance  from  the  wall, 

l  =  Kf  >  (/-//) 
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is  derived  from  equations  (1-7),  (1-8)  and  (1-11),  where  V  is  the 
kinematic  viscosity  and  C  is  a  constant.  This  equation  cannot  give 
the  value  of  the  maximum  velocity  at  the  center  when  we  put  R/ 
because  it  has  been  derived  for  the  flow  near  the  wall.  When  the  de¬ 
viation  from  the  actual  value  is  denoted  by  ^  ,  the  maximum  velocity 

Uo  given  by  the  equation 

-ye=-i-/„-^£  +  C+A  •  (HS) 

Then  we  obtain  the  velocity  defect  from  equations  (1-14)  and  (1-15) 
as  follows: 

^°u  ^  A  )  i or  y/R  « I .  (/~I6) 

(ii)  Near  the  center  of  a  pipe,  or  (/-?/«)«/ 

In  this  case  the  mixing  length  is  represented  approximately  by 
the  following  expression 

•j  — a? -4* (I --%;)*  ■  0-/7) 


From  equations  (1-7),  (1-8)  and  (1-17)  the  differential  equation 
for  the  mean  velocity  results  in  the  form 


dJJ  V/~  y/ft  Uy 
dy  a.1 -4*0- MR?  R 


10 


Using  another  expansion  of  "td/tiAT^X  , 


11 


another  approximate  formula  of  the  velocity  defect  can  be  derived. 


0-2Z) 


The  four  numerical  constants  H  ,  4- ,  C  and  X  have  to  be  determined 
so  that  Equation  (1-16)  for  y/R  «  /  may  be  smoothly  joined  to 
Equation  (1-20)  for  (/-  y/R)  <  <  /  at  some  point. 

By  equating  right-hand  sides  of  these  equations  and  their  de¬ 
rivatives  ,  respectively,  we  have 


and 


X 

ft 


J0£L. 

7 


0-W) 


Now  we  have  four  equations  (1-12),  (1-13),  (1—23)  and  (1-24)  for 
the  four  unknown  constants,  provided  that  the  value  of  X  is  taken 
to  be  0.40  as  usual.  It  is  seen  from  Equation  (1-24)  that  the  value 
of  b/a  changes  rapidly  between  j0£»£3|and  0,  3/  when  X  is 
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substituted  by  Equation  (1-13) .  If  the  point  of  Tf/R  ™  0.  302  J 
is  chosen  as  the  joining  point,  we  obtain  the  values  of  constants  as 
follows: 


and 


<£  a.  0,14  , 

cz  o,o  6  , 


0'2$) 


Putting  these  values  into  equations  (1-10),  (1-16),  (1-20),  (1-21) 
and  (1-22),  the  results  come  as  follows: 

■J-  0M—0.0g(/-£f-O.O6(l-£f'.  o-u) 

=  -S.7S loflJL  JrO.ntS  U-27) 

Ut  K 

for  1/R<0.30 , 

= io. -  ta'jo.stfjhfjj  u-2t) 

for  Y/R  >0.30  , 
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and  approximately 

-^=4.76  (i-£f{i+o.244(l-£f) ,  0-2  f) 

or 

= -S-.43|ln  [/-/■  Hjh%]  +i.HjhJc}-  ( 1-30) 

The  expression  for  mixing  length  is  the  same  as  that  of  Nikuradse.  The 
value  of  A  in  the  equation  for  velocity  distribution  near  the  wall  is 
smaller  than  that  of  Hinze.  The  classical  formula,  which  was  empirically 
deduced  by  Darcy*4, 

J^*=$.08(l-£/Z  ,  0-3/) 

being  the  first  approximation,  Equation  (1-29)  maybe  regarded  as  a 
second  approximation.  Equation  (1-30)  is  akin  to  the  universal  velocity 
distribution  deduced  from  Karman's  similarity  law. 

^ - f  \fh%]  ■  0-32) 

In  Fig.  2,  Equations  (1-27)  and  ( 1—28)  are  plotted  to  be  compared  with 

the  experimental  results  taken  from  the  date  of  channel  flow11  and 

15  16 

pipe  flow  by  Laufer  and  of  pipe  flow  by  Komatsu  .  It  also  contains 

the  velocity  distribution  by  the  Law  of  Wake  deduced  by  Coles,  which 

is  modified  with  the  outer  wall  function  by  Cornish.  The  experimental 

results  seem  to  depend  appreciably  on  Reynolds  Number,  and  there  is 

a  slight  difference  between  the  mean  velocity  in  channel  flow  and  that 
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in  pipe  flow.  The  velocity  distribution  by  Cornish  coincides  fairly 
well  with  the  channel  flow,  while  the  present  distribution  curve  passes 
through  the  midst  of  all  measured  points. 

1-3 .  The  mean  velocity  distribution  in  the  outer  layer  of  a  turbulent 
boundary  layer 

As  mentioned  in  1-1,  if  the  velocity  in  the  turbulent  part  of  the 
intermittently  turbulent  region  is  assumed  to  be  expressed  by  the 
formula  for  a  fully  turbulent  pipe  flow,  Equation  (1-28),  the  mean  velocity 
distribution  in  the  outer  layer  of  a  turbulent  boundary  layer  may  be  cal¬ 
culated  by  equations  (1-2)  and  (1-28).  That  is, 

I0-868  r\MT'{o-uiJh%)  -  .  u-33) 

where  is  the  intermittency  factor  given  by  Equation  (1-4).  It  can 
be  seen  in  Fig.  1  that  the  velocity  distribution  was  not  improved  by 
Equation  (1-33) . 

Now  it  appears  that  the  intermittent-mean-velocity  hypothesis 
is  not  valid  in  the  outer  region  of  turbulent  boundary  layer.  In  the 
following  chapter,  another  hypothesis  will  be  proposed. 


CHAPTER  2 


THE  INTERMITTENT-EDDY- VISCOSITY  HYPOTHESIS 
2-1 .  Introduction 

In  the  preceding  chapter,  the  shearing  stress  was  assumed  to 
be  given  by  Equation  (1-7)  according  to  the  momentum  transfer  theory, 
and  to  be  distributed  linearly  as  shown  in  Equation  (1-8).  The  assump¬ 
tion  is  made  for  a  pipe  flow,  not  for  a  boundary  layer  flow.  Besides, 
the  concept  of  the  intermittent-mean-velocity  hypothesis  is  not  based 
on  any  other  sound  foundation  from  the  dynamical  point  of  view  than 
from  the  phenomenological  one.  Now  we  shall  start  from  the  fundamental 
equations  of  motion  for  the  boundary  layer  flow.  Concerning  the  ex¬ 
pression  of  shearing  stress  in  terms  of  mean  velocity,  there  are  several 
ways  as  found  in  textbooks.  The  first  of  them  is  the  mixing  length 
theory,  including  the  momentum  transfer  theory  and  the  vorticity  trans¬ 
fer  theory.  The  second  is  the  mechanical  similarity  rule  by  Karman. 
Prandtl's  new  assumption  for  the  eddy  viscosity  is  related  to  free  turb¬ 
ulent  flow.  The  constant  eddy  viscosity  concept  by  Boussinesq  comes 
last,  though  it  is  the  oldest  one,  because  it  is  newly  applied  by  Townsend 
and  Hinze.  According  to  the  experimental  results  by  Klebanoff,  the 
eddy  viscosity  in  the  turbulent  boundary  layer  is  not  constant,  but  it 
has  a  maximum  value  near  the  layer  of  “%/&  —0>3  ,  and  it  changes 
slightly  between  y/s  =  0.2  and  0,  S~ ;  whereas  outside  of  the  region 
it  falls  rapidly  to  zero.  The  assumption  of  constant  eddy  viscosity  may 
be  permissible  only  for  the  region  of  y/s<o.s  .  In  the  region 
where  the  eddy  viscosity  changes  rapidly  to  zero,  the  intermittency  of 
flow  pattern  can  be  observed.  It  is  shown  by  Hinze  that  the  eddy  viscosity 
divided  by  the  intermittency  factor  is  nearly  constant  throughout  the  outer 
region,  in  the  same  way  as  it  is  shown  by  Schubauer  that  the  turbulent 
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energy  distribution  divided  by  the  intermittency  factor  shows  a  fairly 
good  agreement  with  in  pipe  flow.  This  factor  shows  a  fairly  good 
agreement  with  in  pipe  flow.  This  fact  seems  to  suggest  the  intermittent- 
eddy-viscosity  hypothesis  proposed  in  the  following  section. 

2-2 .  The  intermlttent-eddy-viscoslty  hypothesis 

The  shearing  stress  in  turbulent  flow  is  given  in  terms  of  the 
eddy  viscosity  as  follows: 

j  ’  (2-t) 

where  £  denotes  the  average  eddy  viscosity  in  a  layer  in  the  outer 
region,  and  not  having  a  constant  value,  but  being  a  function  of  y/S". 
According  to  the  suggestion  mentioned  before,  it  is  assumed  that  the 
eddy  viscosity  has  a  constant  value  £a  in  the  turbulent  region,  while 
it  is  zero  in  the  non- turbulent  region,  and  the  eddy  viscosity  of  £0 
and  of  zero  appears  intermittently  at  a  specified  point.  This  is  the 
intermittent-eddy-viscosity  hypothesis  and  it  is  written  as 

<5  =  to  >  (z-z) 

where  f  is  the  intermittency  factor.  The  same  concept  as  this  is 
suggested  by  Hinze.  The  shearing  stress  in  the  outer  layer  is  now 


(2-3) 


In  a  fully  turbulent  flow,  the  shearing  stress  is  mainly  due  to  the 
turbulent  fluctuation  and  is  expressed  as 
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T - f  u-v  • 


(2-4) 


where  U  and  V"  denote  the  turbulent  velocity  fluctuations  in  X-  and 
y- direction,  respectively.  From  equations  (2-3)  and  (2—4)  we  have 

( i-s ) 

We  can  verify  this  relation  by  making  use  of  the  experimental  results 
of  2 Hr  ,  t  and  u  ,  which  were  measured  independently.  The  turbu¬ 
lent  shear  stress  is  plotted  against  the  velocity  gradient  in  Fig.  3. 
According  to  the  constant  eddy  viscosity  concept,  the  shear  stress 
should  be  linear  to  the  velocity  gradient  itself.  The  relations  between 
the  shear  stress  and  the  velocity  gradient  is  shown  with  small  crosses 
in  Fig.  3.  It  can  be  seen  in  the  figure  that  the  assumption  of  constant 
eddy  viscosity  is  valid  between  =  0-Z.  and  O’  5  /  but 

that  the  assumption  should  not  be  applied  beyond  y/S^O.S.  . 
When  the  shear  stress  is  plotted  against  the  velocity  gradient  multiplied 
by  the  intermittency  factor  iT  ,  instead  of  the  velocity  gradient  itself, 
which  are  shown  with  small  circles  in  the  same  figure,  it  appears  that 
the  linearity  between  them  is  well  confirmed  all  over  the  outer  region. 
Thus  the  intermittent-eddy-viscosity  hypothesis  seems  to  be  valid  in 
the  outer  layer  of  turbulent  boundary  layer  flow. 

In  Fig.  3,  the  shear  stress  is  given  in  non-dimensional  form 
divided  by  the  square  of  friction  velocity,  the  velocity  is  given  in  non- 
dimensional  form  of  velocity  defect  divided  by  friction  velocity,  and 
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the  distance  from  the  surface  is  divided  by  the  boundary  layer  thickness. 
Therefore,  the  slope  of  the  straight  line  gives  the  eddy  viscosity  in 
non-dimensional  form,  that  is,  e,  /s  ut  ,  and  the  value  is  0.062 
for  this  case. 


2-3 .  The  distribution  of  mean  velocity  and  shear  stress  in  the  outer  layer 
of  turbulent  boundary  layer 

The  equation  of  motion  is  given  for  the  mean  velocity  in  turbulent 
boundary  layer  with  zero  pressure  gradient  as  follows: 


1 1  dU  ,  w3Li  _  I  2T 
U3X  v  3^  j>  ay.  ' 


(2-6) 


or  in  terms  of  velocity  defect, 

U.(/-^)  §z(U.-U)  +  V~(u.-U) - j-lf-  (2-7) 


In  the  outer  layer  we  can  assume  that 


*Er«  I 


> 


and  by  the  order-of-magnitude  procedure  we  arrive  at  the  approximate 
equation  of  (2-7) , 


]_ 

f  3?  ' 


(2-8) 


u^(u.-u) 


(2-m 
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where  the  prime  means  differentiation  with  respect  to  £  .  Except 
near  the  edge  of  boundary  layer,  ~f  and  are  of  the  same 

order  of  magnitude.  As  shown  in  the  Appendix, 

_ l_  d.  Ur  ^  /  J_  dS 

UT  dx  S  dx 


Hence,  it  is  permissible  to  neglect  the  first  term  on  the  left-hand 
side  of  Equation  (2-11),  which  then  becomes 


(2 -12) 


Since  the  left-hand  side  of  this  equation  is  a  function  of  t?  alone, 
the  equation  can  be  correct  only  if 


U0S  d.5 _ 

ec  dx  *  ' 


(2-13) 


that  is,  a  constant,  independent  of  X  .  Integrating  Equation  (2-12), 
we  have 


From  equations  (2-9)  and  (2-10) 


(2-14) 


and  since  ^  —  f  and  f  mm  f  when  £ 


[*n 


Ur* 

to 


Hence 


ln(U£) 


-  dlidf  ’ 


(1-1 B) 


=-f'j  <*->« 


Integrating  this  equation  from  infinity  to  ^  ,  we  obtain  the  velocity 


distribution  function  as  follows: 


’  (2~l7) 

since  "f  tends  to  zero  when  £  goes  to  infinity.  For  the  shear  stress 
distribution  we  have  from  equations  (2-14)  and  (2-16) 

]•  (z"l8) 


In  order  to  determine  the  values  of  two  constants,  oi^  and  UtS/e, 
using  Equation  (2-1S) ,  lo^(—  is  plotted  against  fh/rUf 

in  Fig.  4,  when  the  measured  values  of  2T  and  "f  by  Klebanoff  were 


used.  It  can  be  seen  that  the  plotted  points  are  on  a  straight  line. 
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It  suggests  validity  of  the  relation  given  by  (2 
two  constants  obtained  are: 

,  or  £ 

and 

d1  =  4-so  ■ 

This  value  for  non-dimensional  eddy  viscosity,  U/Ur  s  ,  is  in 
good  agreement  with  that  obtained  before  from  the  shear  stress  in 
Fig.  3.  The  mean  velocity  distribution  of  Equation  (2-17)  with  the 
constants  of  (2-19)  is  shown  in  Fig.  1.  The  calculated  values  show 
an  excellent  agreement  with  those  measured  by  Klebanoff.  The  velocity 
distribution  by  the  law  of  wake  modified  with  the  outer  wall  function 
by  Cornish  is  also  shown  in  the  same  figure.  This  distribution  will 
be  discussed  in  the  next  section. 

The  shear  stress  distribution  calculated  by  Equation  (2-18)  with 
the  value  of  dL*'  of  (2-19)  is  shown  in  Fig.  5,  and  it  is  also  in  good 
agreement  with  that  measured  by  Klebanoff. 

2-4 .  The  wake  function  and  the  outer  wall  function 

Coles2  has  proposed  that  the  mean  velocity  distribution  may  be 
written  in  the  form 


-15).  Values  of  the 


>0-062 


(2-W) 


(2-20) 


23 


which  is  now  known  as  the  law  of  wake,  and  ifriy/S)  is  called  the 
wake  function.  TT(X)  denotes  the  profile  parameter.  H  is  an 
empirical  constant  in  the  logarithmic  law ,  that  is 

f(JtMr)  _  JL  +  C  ,  (2-21) 


where  C?  is  another  empirical  constant.  Coles  used  the  numerical 
values  of  ) (=0.40  and  C  — S-IO .  The  profile  parameter  TT  (x)  is 
given  the  value  of  0.55  for  a  flow  at  constant  pressure.  Then  the 
velocity  profile  is  written  as 


J=L  =  J—  lnl^-r  +  5. io  +  £££  mr(JL) , 
UT  040  1/  +  1 


or  in  the  form  of  velocity  defect 


+  ffH2— 


since  the  wake  function  is  normalized  so  that  it  has  the  value  $t( l)*"Z. 
From  this  equation  we  have 


1 -  ■ 2 -  Hr +  sfe  W)  ■  (2-Z2) 


By  using  Equation  (2-17)  for  the  velocity  defect,  the  value  of  w(y/S) 
can  be  calculated.  The  result  is  shown  in  Fig.  6  with  that  of  Coles. 
The  wake  function  by  Coles  is  seen  to  be  smaller  than  the  present 
result.  Besides,  the  calculated  wake  function  has  the  maximum  at 
about  y/$  as  O.^  ;  whereas  Coles’  wake  function  has  no  maximum 
value.  This  leads  to  the  non-zero  slope  at  the  edge  of  boundary 
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layer,  as  was  pointed  out  by  Cornish. 

The  outer  wall  function  has  been  proposed  by  Cornish3  to  im¬ 
prove  Coles'  description  of  the  velocity  profile,  that  is, 

■u  ” T1"' +  C  +  -F «*#>  +  “>(#)  ’  (2-23> 

where  Qj(y/$) is  the  outer  wall  function.  The  velocity  defect  is  written 
as 

since  U)(l)=0 •  With  X^O.40  and  TT  =  Q.SS  #  this  equation 
becomes 

from  which  we  have 

"«#>— /n  f  +  ff  C2~^J  -  TT  •  (2'24) 

By  using  the  values  of  Coles  for  the  wake  function  and  Equation  (2-17) 
for  the  velocity  defect,  the  outer  wall  function  can  be  calculated, 
and  it  is  shown  in  Fig.  7.  The  difference  between  the  outer  wall 
function  of  Cornish  and  the  present  result  is  very  small  except  between 
y/s=o,e  and  0.8.  It  should  be  noted,  however,  that  the  values  of 
constants  used  by  Cornish  were  K 0.41  Z  and  ZV/K=2.5  . 
so  that  the  value  of  TT  was  0.515  .  With  these  values  the  discrepancy 
will  be  seen  to  become  larger. 
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2-5.  The  shape  factor 

The  shape  factor  is  defined  as  the  ratio  of  displacement  thick¬ 
ness  and  momentum  thickness  of  boundary  layer.  Velocity  defect 
is  similar  if  it  is  divided  by  friction  velocity,  not  by  free  stream 
velocity.  Accordingly,  the  shape  factor  is  a  function  of  friction 
velocity.  The  relation  between  them  is  written  as  follows: 


H-4-- 


«  /- 


fcu r 

1,  Up 


where 


(2  -2S) 


(2-26; 


lfU0  ’ 


8  r  Ur  _T  Ur 

S  L'Ue  12  U* 


order  to  estimate  the  values  of  X/  and  I2  >  we  use  the 


logarithmic  law 


•^uT  =  —  A  In  ?  +B 


(2  -27) 


with  A=Z5  and  B— 2*3  from  ^  =  Q  to  and  the  present 

formula  of  equation  (2-17)  with  (2-19)  from  ^=*0,/  and  /  , 

and  we  obtainI/=3.676’  ,  1^=24,50/  andJa  /li=6.G7-  with 
this  value  of  Xa/r, .  H  is  plotted  against  Ux/Uo  In  Fig .  8 . 
This  result  is  seen  to  be  in  a  better  agreement  with  the  measured 
values  than  that  of  Hama  with  the  value  of  Xa/jQ  *  /  .  It  may 


be  noted  that  if  the  logarithmic  formula  (2-27)  is  used  throughout 
the  boundary  layer,  the  ratio  is 

Ia_  2AV2AB  4-B* 

1/  A  +  & 

and  the  value  Il/I I  -  6.  I  is  obtained  with  A  =2, 5  and  8=2.3. 
If  we  use  Equation  (2-17)  only,  we  obtain  Tjl,-S.S  7  .  With 
Cornish's  outer  wall  function,  almost  the  same  result  as  mentioned 
above  is  obtained  as  shown  in  Fig.  8. 


2-6.  Eddy  viscosity  and  mixing  length 

Clauser  noticed  that  the  constant  pressure  turbulent  profile 
dropped  so  abruptly  at  the  wall  as  to  appear  to  extrapolate  to  a  non¬ 
zero  velocity  at  the  wall,  and  that  this  characteristic  shape  of  the 
turbulent  profile  came  from  the  circumstances  that  the  laminar  sub¬ 
layer  next  to  the  wall  and  the  flow  adjacent  to  it  had  a  lower  viscosity 
than  the  eddy  viscosity  prevailing  in  the  main  body  of  the  turbulent 
flow.  He  attempted  to  simulate  a  turbulent  boundary  layer  profile 
to  a  laminar  profile  with  a  slip  velocity  on  the  wall  and  with  an 
appropriate  eddy  viscosity,  instead  of  a  laminar  viscosity. 

If  the  present  formula  of  Equation  (2-17)  is  extrapolated  to  the 
wall,  the  velocity  at  the  wall  or  the  slip  velocity  Us  Is  obtained. 


Uo-Ut 

U* 


nr  UxS 

2 i1  £. 


(2-2® 


which  is  calculated  to  have  a  constant  value  of  9.54  with  (2-19). 
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As  shown  in  preceding  sections ,  eddy  viscosity  cannot  be 
regarded  as'constant  in  a  turbulent  boundary  layer,  but  it  is  nearly 
constant  in  the  turbulent  part  in  the  outer  layer.  The  eddy  viscosity 
in  the  turbulent  part  is  written  from  Equation  (2-5)  as  follows: 


and  6.0  /us  is  0.062  as  mentioned  before.  From  Uf»1.85  ft/sec. 
and  S»3.0  in.  of  Klebanoff's  data,  we  have £,=  0.0287  ftVsec. 

Since  ,  from  Equation  (2—26)  and  Iysss3.675, 

we  have  another  expression  of  £0  , 

£0=  0-0168  US*  (2-21) 

From  this ,  we  have  the  eddy  Reynolds  number 

R  =^=S1-S  (2-3 0) 

^  Co 

The  eddy  viscosity  by  Claus er  is 

£  —  rfU.  S* 


and  0(  is  universally  0.018.  Therefore,  the  eddy  Reynolds  number 
is  U.S*/£~S6  .  It  is  noted  that  the  present  results  coincide 
rather  remarkably  with  those  of  Clauser  in  spite  of  the  difference  in 
definition  of  £®  and  £  . 
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In  definition  of  given  by  Equation  (2-3),  the  transfer  of 
mean  flow  momentum  is  assumed  implicitly  to  be  due  to  the  small 
eddy  or  gradient  type  of  transfer,  the  large  eddy  of  convective  type 
of  transfer  being  neglected.  From  the  constancy  of  shown  before, 
it  may  be  assumed  that  the  turbulent  part  in  the  outer  layer  is  com¬ 
posed  of  small  eddies  of  the  same  size  and  of  the  same  strength. 

In  other  words ,  the  motion  of  a  fluid  element  may  be  regarded  as 
random  walk  of  flight  velocity  V"  and  of  flight  length  G.  .  Then  the 
eddy  viscosity  is  written  as  follows: 


where  the  flight  length  H  corresponds  to  the  mixing  length.  Accord¬ 
ing  to  the  measurement  by  Klebanoff,  the  y-component  of  turbulent 
intensity  ■y-/tj^iivided  by  the  intermittency  factor  has  roughly  a  con¬ 
stant  value  between  0.035  and  0.043.  If  we  take  the  value  vAj„=0M, 
we  have  from  Equation  (2-29) 

J~  =  0-84- 
8* 

With  Equation  (2-26)  and  the  value  ofI/=»  3.765  and Ur/l-t3* ® •  03 7 , 
the  mixing  length  is  obtained  in  terms  of  the  boundary  layer  thickness, 
that  is 


JL  =  0-II7 
S* 


(2  -31) 
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It  can  be  seen  that  the  mixing  length  is  small  compared  with  the 
boundary  layer  thickness .  This  result  is  consistent  with  the  assump¬ 
tion  of  gradient  type  of  transfer  of  the  mean  flow  momentum.  The 
effect  of  large  eddies  is  taken  into  account  through  the  intermittency 
factor. 


CONCLUSION 

A  new  formula  for  the  mean  velocity  distribution  of  a  fully 
turbulent  pipe  flow,  which  was  derived  with  reference  to  the  mixing 
length  measured  by  Nikuradse,  agrees  well  with  experimental  results. 
The  mean  velocity  profile  in  the  outer  layer  of  a  turbulent  boundary' 
layer  was  calculated  by  the  intermittent-mean-velocity  hypothesis, 
using  the  new  formula  for  pipe  flow,  and  the  result  did  not  agree 
with  the  experiments . 

It  has  been  shown  by  study  of  measured  results  by  Klebanoff 
that  the  eddy  viscosity  in  turbulent  region  of  outer  layer  is  nearly 
constant  and  that  in  non-turbulent  region  is  zero,  so  that  the  turb¬ 
ulent  shear  stress  is  expressed  in  terms  of  the  average  eddy  viscosity 
which  is  the  product  of  the  constant  eddy  viscosity  and  the  in*.  - 
mittency  factor.  It  is  suggested  to  call  this  concept  the  intermittent- 
eddy-viscosity  hypothesis.  The  velocity  profile  and  the  shear  stress 
distribution  in  the  outer  layer  deduced  from  this  concept  are  in  good 
agreement  with  the  experimental  results . 

Coles'  wake  function  should  be  modified,  and  for  this  purpose 
Cornish's  outer  wall  function  gives  a  little  larger  values  than  those 
calculated,  though  the  difference  is  very  small. 

The  relation  between  the  shape  factor  and  the  friction  velocity 
is  well  described  with  the  present  velocity  profile. 

The  eddy  Reynolds  number,  UoSyeo  ,  has  the  value  of  60, 
and  the  mixing  length  has  a  constant  value  of  order  of  one-tenth  of 
boundary  layer  thickness,  so  that  the  gradient  type  transfer  of  mean 
flow  momentum  can  be  assumed. 
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APPENDIX 


Assuming  that  the  velocity  is  a  function  of  Ut  y/v, 


and  substituting  S  for  y  ,  we  obtain  the  freestream  velocity  Uo, 

-%=F(^)  (2) 

Differentiation  of  Equation  (2)  with  respect  to  X  yields 


where  p^enotes  the  derivative  of  F  with  respect  to  U^S/i/  .  Since 

2 

the  local  skin  friction  coefficient  Cx*m2.(UT/U0)  ,  Uo/Ux  becomes 
very  large  at  large  Reynolds  numbers.  If  the  velocity  is  assumed  to 
have  the  logarithmic  form 

ft-Fffi-Aln**  +  6  (+> 

we  have  the  following  relation  by  differentiation 

1 _ £_  =  a 

F'  m 


(S) 


Strictly  speaking,  this  relation  cannot  be  valid  in  the  outer  layer. 
However,  the  order-of<-magnitude  relation  will  remain  unchanged. 
Therefore,  S  ,  may  be  of  the  order  of  unity,  because  A 

has  the  value  of  about  2.5.  Accordingly,  the  second  term  in  the 
bracket  on  the  right-hand  side  of  Equation  (3)  can  be  neglected.  In 
other  words,  it  is  very  small  in  comparison  with  the  left-hand  side. 


Uo  d.  (  Ur\  /  ✓  _]_  dLS 
Uxdz'Uo'  ^  S  S  dx: 


(6) 
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Figure  8 


